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: ABSTRACT 
^ . 

^ ■ We argue that primordial black holes in the early universe can provide an efficient resolution 

Q^ ■ of the Brustein-Steinhardt moduli overshoot problem in string cosmology. When the universe 

is created near the Planck scale, all the available states in the theory are excited by strong 
O '• interactions and cosmological particle production. The heavy states are described in the 

low energy theory as a gas of electrically and magnetically charged black holes. This gas 
of black holes quickly captures the moduli which appear in the relation between black hole 
masses and charges, and slows them down with their vevs typically close to the Planck scale. 
^ . From there, the modulus may slowly roll into a valley with a positive vacuum energy, where 

inflation may begin. The black hole gas will redshift away in the course of cosmic expansion, 
^ : as inflation evicts the black holes out of the horizon. 
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Inflation [1] is at present the best framework for explaining the origin of our universe. 
However there stifl are technical difficulties with implementing inflation in fundamental the- 
ory. Inflation needs light scalar fields with very fiat potentials and positive vacuum energy, 
which are hard to obtain from first principles. Even if such degrees of freedom are present, 
there is still the problem of arranging favorable initial conditions [2] , which ensure that the 
vacuum energy controlled by the light scalars dominates over other energy sources in the in- 
fiationary patch. One approach towards resolving this problem is the idea of eternal inflation 
[3, 4, 5] (for a recent review see [6]), which posits that after inflation starts in some regions 
of a huge "metauniverse" where the environmental conditions favor it, it will produce many 
universes as big as ours. 

Recently, there has been progress in the direction of implementing eternal inflation in 
string theory. Backgrounds where eternal inflation can occur have been constructed in flux 
compactiflcations [7, 8, 9]. These solutions naturally belong to the landscape of string vacua 
[10, 11, 12, 13]. In fact, many of the older cosmological solutions in supergravity limits of 
string theory without future singularities [14] could be fltted at the foothills of the landscape, 
close to the supersymmetric limits. The landscape picture provides for many valleys where 
inflation can happen, with some of the moduli flelds as inflatons. This would realize the old 
hope that string moduli may be inflatons [15]. However, this approach may also resurrect the 
moduli overshooting problem, discussed by Brustein and Steinhardt [16]. Depending how 
they start, the moduli flelds may acquire a lot of kinetic energy on the approach towards the 
inflationary minima and overshoot them, running off towards the supersymmetric regions of 
weak coupling. There the spacetime decompactifles and new dimensions of space open up. 
Thus in order to flnd a phenomenologically viable cosmology, it is essential to flnd a way to 
gently lower the moduli to where inflation can occur and the world remains four-dimensional 
for a long time. 

In this framework, the universe may be described as a compact space, with flat or nega- 
tively curved spatial slices [17] (see also [18]). It is created close to the Planck scale. Three of 
the spatial dimensions grew very big during inflation. The modulus which plays the role of 
the inflaton could start in the inflationary valley of the potential like in topological inflation 
[9, 17]. However in some regions modulus may have started high up the slope, and normally 
it would have overshot the minimum as it rolled down; if instead there is enough radiation, 
the modulus could still get captured in the minimum thanks to the extra cosmological fric- 
tion [19]. The main ingredient of this mechanism has been observed in string cosmology 
earlier in [20, 21]. In the backgrounds dominated by conformal matter where T^^ = the 
source terms arising from the gravitational couplings of the moduli to matter vanish and 
moduli are quickly kinetically damped [22, 23]. Their energy density dilutes as ph ^ for 
the homogeneous mode or as pj ~ 1/a^ for inhomogeneities [24], where a is the cosmological 
scale factor. Some aspects of more phenomenologically motivated, but similar, mechanisms 
have also been discussed in [25, 26, 27]. 

We will present here a different mechanism for moduli entrapment. The idea is to use 
the source term coming from the interactions of the moduli with very heavy, nonrelativistic, 
charged states in the early universe to enhance the attraction of the inflationary basin. To 
illustrate the mechanism consider the following analogy with a marble and a bathtub. If 
one drops a marble alongside the steep wall of the bathtub, the marble will have acquired 
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so much kinetic energy that it will fly over any of the small dimples in the tub's floor and 
continue rolling around for a long time. However, if the tub is first filled with water, which 
slowly drains away, then the marble will loose its kinetic energy very rapidly thanks to the 
interactions with the water molecules. It will be captured quickly in a dimple on the bottom 
near the place of original impact. Then water will drain away, analogously to the redshifting 
of the cosmological matter contents, leaving the marble, or the modulus potential, in control 
of the evolution of the universe. We will show that a gas of heavy states in string theory, 
described in the low energy supergravity limit as extremal black holes, can enforce the same 
effect on the runaway moduli of string compactifications. In this way, the gas of black 
holes, with high initial energy density, can provide a new mechanism for overcoming the 
Brustein-Steinhardt overshoot problem. 

In what follows we will work with a single dilaton-like scalar modulus, assuming that 
the theory has been compactified to four dimensions already. We believe that the same 
mechanism can be realized in more complicated scenarios. Then in the approach of [17], 
the universe emerges somewhere on the landscape with an energy density close to the string 
scale. We will consider what happens with the regions where the modulus is not immediately 
placed in the inflationary plateau. Instead it starts far from the minimum along a steep slope 
and acquires a large kinetic energy. Because the spatial curvature of the universe is either 
zero or negative, it does not immediately collapse, but expands under the influence of the 
dominant sources of stress-energy. This separates the flatness problem from the Planck 
scale physics, and allows for the possibility that inflation may begin well below the Planck 
scale. At very high energies where the universe is born, strong interactions, and specifically 
gravitational particle production [28, 29] will excite all the modes in the spectrum of the 
theory. This includes the heavy states in the theory, with masses above the string scale. 
When their masses exceed the Planck scale, such states, with some given quantum numbers, 
are described as BPS black holes in the supergravity hmit. Such states are known to locally 
fix the scalar field vevs on the horizons to values completely independent of the asymptotic 
geometry [30]. 

The masses of these states depend explicitly on the expectation value of the scalar at 
infinity. To see this, recall the example of the BPS black hole solutions in heterotic string 
theory, described by the bosonic sector of the effective supergravity action [31, 32, 33] 

s^l d'x^g ^ (it: + {d^f - \fI:) . (1) 

A generic BPS soliton of (1) is described by an extremal black hole configuration with mass 
M., electric charge Q and magnetic charge which are related by [31, 32, 33] 

M = l(|Q|e^/2 + |7^|e-^/2), (2) 

where (p is the dilaton vev infinitely far away from the hole, determining the asymptotic value 
of the string coupling gs = exp{ip/2). This parameter is left completely undetermined by 
the local black hole dynamics. In the supersymmetric limit it is a modulus which can take 
any value because it is a flat direction of the theory. The charges Q and V are quantized in 
the units of the string scale ~ 1/is- Similar situation persists for the case of other moduli 
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fields. If, for the sake of simplicity, we assume that all but one of the moduli are stabilized 
in the compactification to four dimensions (such as in [7, 35]), the effective 4D action which 
describes the light modes of the theory is, in the Einstein frame, [31] 

S = J d'x^g - \{d^f - Y^^'-^Fl) , (3) 

where g2 is the modulus coupling, and k — l/Mpi is the inverse Planck mass in 4D. A generic 
BPS soliton of the theory (1) is described by an extremal black hole configuration with mass 
M., electric charge Q and magnetic charge V. In general, the solutions for the arbitrary 
charge assignment for Q, V are not known exactly for an arbitrary value of g2 [31]. However, 
it is straightforward to solve (3) for the black hole solutions with only electric, or only 
magnetic charge, in the extremal limit. They are given by zero entropy, zero temperature 
configurations^ of the same global structure as those with g2 = -\/2 [31]. Their masses and 
charges are related according to^ 

K K 

where again (p is the scalar vev infinitely far away from the hole, determining the asymptotic 
value of the coupling constant g = exp(g2fi;0/2), and can take any value in the supersymmet- 
ric limit. In what follows we will work with the models where the Planck scale and the string 
scale are close together, so that asymptotically g ~ because typically the minima of 

the nonperturbative scalar modulus potential appear in this regime. Appropriately, as we 
will see in what follows, the black hole toy model works as a moduli stopper most efficiently 
precisely in this regime. 

The scalar minima in the mass function (2) have prompted a suggestion [36] that moduli 
could acquire effective potentials generated by the virtual black hole loops, which can sta- 
bilize them. This idea is very interesting because it may incorporate some nonperturbative 
quantum gravity effects, even though is difficult to reliably compute them. We still don't 
know quantum gravity Feynman diagrams for black hole states. Instead of following this 
route, we will focus on a different application of the minima of (2). We will consider an early 
universe environment where there is a large number of black holes with scalar-dependent 
masses as in (4). Once there is a net energy density of them, with both electric and mag- 
netic charges, their interactions with the scalar will produce a damping effect, trapping the 
scalar to vevs of order of the Planck scale (after we normalize the scalar field canonically) . 
The electric charges block the modulus from running off to the strong coupling, where they 
are heavy, and the magnetic charges block it from going to the weak coupling, where they 
become heavy. Thus as long as the charges hang around the modulus cannot easily escape to 
either side. This is akin to the environmentally induced stabilization effects of [25] . In what 
follows we will only consider the effects of the black hole gas populated with only electric 
and only magnetic charges. The dyonic solutions should also exist, and contribute as well. 

"'^Thcy arc null naked singularities, unlike the Schwarzschild solutions whose interplay with scalars has 
been studied in [34]. Hence they will not Hawking-evaporate. This is why they are identified with heavy 
states of the theory. 

^We have absorbed any irrelevant numerical factors into the definition of couplings. 
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However in the regime where we will work, dyons will be generally heavier than monopoles. 
Because we assume a nearly thermal initial population of heavy states, we expect that the 
dyon contributions and effects will be exponentially suppressed relative to the monopole 
states, and we can safely ignore them. 

How do we describe this ensemble of black holes? It is well known that the extremal black 
hole states with like charges can be arranged in static arrays because the electromagnetic 
repulsions can neutralize the gravitational attraction between them. In the early universe 
these charges would neither remain static, nor would they all be of the same sign. Gravity 
respects gauge symmetries and so charge conservation would require the net charge to be 
zero. In general, the black hole states will be initially produced by strong (gravitational) 
interactions near the string scale [28, 29] both as "particles" and "antiparticles" , with oppo- 
site charges, carrying both electric and magnetic charges. Hence we will approximate their 
initial abundances by the nearly thermal distribution law, which should be a sensible order 
of magnitude estimate. Systems of such black holes would be dynamical, since the forces 
between black holes would not cancel exactly, due to the presence of charges with both signs. 
Furthermore, outside of the black holes the space would not be in the vacuum, but in a some, 
roughly thermally, excited state. It is clear by CPT invariance (4) and equivalence principle 
that both "particles" and "antiparticles" contribute equally to the stress-energy. States from 
all mass scales will contribute to the total stress-energy. If they meet and annihilate, their 
energy would be released as relativistic particles, which we will ignore later on. The effect 
of radiation on cosmological moduli dynamics has already been considered in [20, 21, 19]. 

While overall unstable, such arrays of black holes may become sufficiently separated so 
that the cosmic expansion prevents their complete annihilation. In this case, away from the 
static limit the leading order electromagnetic interactions can be neglected, and the leading 
order gravitational effects from the charges can be approximated by the stress-energy tensor 
of the fluid of massive particles [37] . Thus since they are heavy states, almost immediately 
after they are produced, surviving black holes will fall out of thermal equilibrium. Therefore 
their leading contribution to the energy density will come from their rest masses. With this 
in mind wc can use the dilute gas approximation for the description of the evolution of these 
particles. This should be reasonable below (but close to) the string scale, and will rapidly 
become better as the universe expands. The initial energy density of black holes should 
be high, but still below the Planck scale, to guarantee the validity of the (super)gravity 
description where we can ignore higher derivative corrections in the effective action. 

We expect that while our approach is an approximation, it should be a reasonable one in 
most of the space when the initial black hole density is not too large. Namely, it is clear that 
close to any individual black hole the averaging procedure which we will embrace below must 
break down, because of the strong fields close to a hole [30]- [36]. In that region, the solution 
is completely controlled by the local dynamics which fixes the modulus to a value specified 
by the charges and renders it insensitive to its value at infinity [30, 36]. However we are 
interested in showing how an initial population of black holes can help resolve the overshoot 
problem [16] of string cosmology, which is the statement that the asymptotic value of the 
scalar is really a zero mode that may run off to the weak coupling regime. If the black hole 
density is not too large initially, so that their individual separation is e.g. >, (9(10) Planck 
lengths (and initially 1/Hq > ipi), then the value of the modulus in the space surrounding 
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them will rapidly approach its asymptotic value. Our approximation should be a reasonable 
approximation controlling the dynamics of the modulus in the interstitial space between the 
holes. Once inflation starts, this approximation will rapidly get better and better as time 
goes on. In a more realistic situation, however, the modulus will have a distribution of values 
determined by the individual black holes. There will be pockets with all kinds of values of 
couphngs locally fixed by hole charges, and those which can be later stabihzed by different 
minima of the landscape potential where infiation can occur will give rise to sibling universes 
with different low energy couplings^. 

We can estimate the initial energy density of black holes as follows. The number density 
of the black hole states which are produced is Boltzmann-suppressed, um oc exp(— AI/Tq), 
where Tq ~ y/HoMpi is the temperature of the universe when it is formed. We will take it 
to be To ^ 1/^5, so that we can reliably use the supergravity limit for the description of 
cosmology. The total energy density which these massive states contribute will be [22, 20] 



where we estimate Ai by (4) and rij^, as indicated above, roughly by the Boltzmann distri- 
bution when A4 > 1/is- Thus the main contribution will come from the lightest black holes. 
From (4) and Tq ~ l/fg we see that the Boltzmann suppression factor of the initial number 
density scales as um oc exp(^— | QlMp^^^^fo) and tim oc exp(^—\V\Mpi£s/goj- For a fixed ini- 
tial value of the coupling constant go, the masses obey M ~ goQMpi ~ VMpi/go 3> l/is, 
and since Mpi ~ l/{gois) they indeed are black holes, M. ^ Mpi. Thus when go is of the 
order of unity, the black holes can be described as a nonrelativistic fiuid, whose pressure 
can be neglected, and whose energy density in an approximately FRW background can be 
estimated from (5) with the help of a saddle point approximation. In fact because the mass 
in (2) is a linear combination of the terms proportional to g and 1/g, assuming a roughly 
thermal initial distribution no ~ e-J^{'t>o)(-s j can s^j^ ^j^g magnetic and electric 
contributions to pbh separately. This yields the initial black hole energy density 



where no is the number density of the lightest black holes in the ensemble, and {q) = 
J2q \ Q\ ^^v{~\Q\Mpiisgo) /no and (p) are ensemble-averaged values of \Q\, \V\ respectively, 
obeying roughly {q) ~ (p) ~ 0{1) when initial value of g is not too far from unity. As the 
universe evolves undergoing cosmic expansion, the energy density of the black hole gas will 
change in two ways: i) it will redshift away according to the usual ~ 1/a^ law for massive 
particles, since we are ignoring the black hole interactions (we expect that this is justified 
for heavy black holes when they are sufficiently separated [37]) and ii) the evolution of the 
modulus will change the coupling and so the mass of the black holes, as is clear from (4). 
When the evolution is smooth, a good approximation accounting for these two effects is to 
represent the energy density as a function of the scale factor and the modulus vev as [22, 23] 



Pbh ^Y.'^mM, 
M 



(5) 




(6) 




a 



(7) 



^We thank Andrei Linde for interesting discussions of this issue. 
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It is convenient to define 0* by exp(g2«;0*) = Then the formula for the energy 

density of black holes becomes 



Pbh 



— cosh( 



(8) 



where Po = y (p) (?) noa^Mpi. Here we allow for the change in pbh via a subsequent adiabatic 
evolution of away from its initial value. In the limits go ^ 1 and go ^ 1 the number 
densities of black holes with large electric charge and with large magnetic charge, respectively, 
are exponentially suppressed. We will comment on these limits later on. 

The energy density of black holes (8) will appear as a source in the Friedmann equation 
governing the evolution of the scale factor of the universe a. In the Einstein frame, this 
equation is 

3Ml,H^ + 3^ = ^ + V{cl)) + PBH + Prad , (9) 

where we have included for completeness the contribution from all rclativistic particles in 
the universe Prad- We stress that we restrict the spatial curvature to be 0,-1 only thanks to 
the arguments of [17], which strongly favor these values at high densities. Thus the collapse 
is averted. The kinetic and potential contributions from the scalar are included as 0^/2 and 
y(0) respectively. Now, to find how the black hole gas affects the modulus, we need the 
equation of motion for its zero mode. A simple way to obtain it is to use the second order 
Einstein equation for the scale factor, and resort to the Bianchi identity to find [22]. The 
equation for a is, at the level of our approximations where pbh — ^ and Prad = Prad/'^i 



a 
a 



6M. 



1 / • 4 \ 

-2 (202 _ 2V{(P) + Pbh + -Prad) , 



(10) 



and hence taking the first derivative of (9), using (8) and ehminating a from (10) yields 



+ 3H(f) + ^ Pbh tanh( 







) = o. 



(11) 



This is the master equation which controls cosmological dynamics of the scalar modulus in 
this problem. 

Let us now turn to the dynamical evolution of the universe with a modulus as governed 
by equations (9), (10) and (11). First ignore the matter contributions. A typical potential 
for the moduli may have some local minima, but is very steep on one side, and asymptotes 
as an exponential function to zero on the other side of the region where the minima lay. An 
example is provided by a potential given in [7] 



V{cP) 



2(72 



+ Wo 



D 



:i2) 



where a — exp (-^2/3 k 0^ . When we consider detailed examples of the scalar evolution 
below, we will take Wo = -10~^, A = 1, a = 0.1, = 3 x 10"^ in units of Mpi following 
[7], that gives a local minimum for the potential at ~ (lO^^Mp^)^. We plot this potential 
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in Fig. (1). The direction of ^ oo corresponds to the weak couphng hmit, where the 
extra dimensions open up and the moduh are the zero modes from the higher-dimensional 
graviton multiplet. If the modulus is trapped in a minimum, the solution will approximate 
a 4D world for a long time, controlled by the tunnelling rate through the barrier, and allow 
inflationary regime [7]. 
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Figure 1: Potential V {(/)) given by eq. (12), in the units of 10~^^Mp;. In the niininiuni V ~ 

However, for a generic initial value of the modulus in the strong coupling, the modulus 
will fall down the potential and the evolution will be dominated by its very large initial 
kinetic energy. The potential will remain subleading because it is too steep: its initial role is 
to just start the phase of kinetic domination. By the time the modulus reaches a minimum, 
it will have acquired too much kinetic energy from rolling down such a steep potential. As a 
result the modulus will skip over the barrier and continue rolling towards the weak coupling. 
This is the origin of the overshoot problem [16]. If such compactifications are to avoid special 
initial condition, they need a dynamical mechanism which will arrest the modulus before the 
flyover. Otherwise, it will be extremely hard to see how to ever keep the modulus around 
the minimum and start early inflation. 

The radiation terms oc prad in (9), (10) can help to slow down the modulus. Their damp- 
ing effects have been studied in [19] and have been observed earlier in [20, 21]. Because of 
the conformal symmetry, the stress-energy tensor of radiation vanishes, and so the radiation 
does not source the modulus. This is manifest in (11) where the radiation terms are absent. 
Then because radiation scales only as p ~ 1/a^ it will overtake the scalar kinetic energy, 
which will damp out as 1/a^, and the modulus will stop. 

We note that the curvature term k/a"^ will work in a similar vein, but even more efficiently 
than radiation, if = — 1. In that case, the curvature term will begin to dominate the 
evolution of the universe soon after it came into existence, and lead to the linear expansion 
a ~ t. This will dilute the kinetic term of the scalar even faster than radiation. Thus if 
either of these terms dominates early on, and the scalar does not overshoot the minimum 
before it stops to control H, it can get caught and inflation could begin. However, in the 
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case of radiation, since the black hole density redshifts only as it will overtake radiation 
within few Hubble times after the beginning, and if initial k is small or zero, the black hole 
gas will play the key role in the evolution of the modulus. In what follows we will therefore 
ignore the Prad and k/a^ terms, and instead focus on the black hole contributions oc Pbh- 

To leading order, the effects of the interactions of the modulus with the black hole gas 
Pbh can be described as an environment-induced mass term, similar to the interactions 
discussed in [25], and more recently to the chameleon fields [38] (sec also [39]). The presence 
of this new mass term oc pbh suggests a new method for capturing the scalar in a potential 
minimum. The interactions with the black hole gas pull it towards 0* where it awaits for 
the contributions from the potential V to overtake the black hole terms. When 0* is divided 
from the supersymmetric region by a barrier in V, the scalar will eventually settle into a 
minimum of V where inflation can take place. It will not overshoot the barrier since its 
kinetic energy will be spent against the black hole mass pumping. However this effect in 
general is not always enough to stabilize the scalar. Indeed, if the gas gets diluted too fast, 
the interactions will become too weak to dissipate the scalar kinetic energy. This can be seen 
immediately from the comparison with a simple particle mechanics with a time-dependent 
mass. For example, if one considers a pendulum with a mass term = a/t'^, one finds that 
the motion is very unstable, and there are modes which move the pendulum away from the 
minimum at zero. However, in our case the black holes do not dilute dangerously fast! 

To see this, consider a borough of the landscape where the modulus (f) measures the 
size of the 6D space of a warped compactification of the type IIB theory according to the 
prescription of [7, 35]. Then the gauge field in (3) can arise from reducing a lOD 3-form such 
that one index is in the internal space, and so the corresponding choice for g2 is g2 = —^2/3. 
Recall that the weak coupling is recovered in the limit — * oo. Consider (11) with these 
parameters when the scalar is close to the effective minimum at 0*. The small perturbations 
'd = (j) — (j)^ obey the linearized version of (11). After black holes start to dominate in (9), 
the solution behaves as a = t^^^ and pbh — ^MpH"^ — AMp/St'^. In this limit the linearized 
equation becomes 

+ 2M + -1? = . (13) 
9 

The solutions for then are 

'^ = '^* + ^ + ^' (14) 

and so thanks to Hubble damping the evolution is stable under small perturbations. This 
is the key ingredient of our mechanism, which guarantees that the scalar entrapment to 
0* = Mpi\ii{gl{p) / {q)) is gentle enough to ensure the loss of the scalar kinetic energy. 

We have studied the capture of the modulus by the black hole gas numerically in order to 
verify the stability of the entrapment in the nonlinear regime. In a typical case, the evolution 
is depicted in Fig. (2). We have taken the initial density of black holes to be ~ 10~^Mp^, a 
number close to, but safely below the Planck scale, in accordance with our assumptions that 
Qq ~ 0(1), the validity of the (super)gravity approximation and that the universe emerged 
near the Planck scale. Since the initial effects of the potential are merely to induce a large 
kinetic energy of the modulus, early on we can replace it by picking a large initial value of 
(j) ~ and set V — Further because of the Boltzmann suppression in (5) for massive 
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states, it is natural to take the initial value of the black hole energy density to be somewhat 
lower than Mp^. In this case the modulus kinetic energy will start dominating the universe. 
However it dilutes fast, and the black hole density starts to contribute quickly. We have 
assumed here that the initial value of g is of the order of unity, so that the location of the 
attractor for is near Mpi. Clearly, the black holes trap the modulus at a value (p^, and the 
modulus kinetic energy quickly becomes subdominant. 




1 102 104 106 108 10^0 10^2 1 102 104 10^ 10^ 10^° 10^2 
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Figure 2: Numerical results for V = 0. Left: evolution of the field (p, in the units of Mpi, as a 
function of time. We depict time in the units of Mpi/ ^ pbh o throughout, and we set 0* = 2Mpi. 
The initial conditions are 4>o = Mpi, (j)o = -Mil, PBH = 10~^ at t = 0. Right: evolution of 
Pbh and the scalar kinetic energy, in the units of pbh o- 



It is possible to verify this behavior analytically. An approximate solution of eqs. 
(9), (10), (11) initially is given by the FRW cosmology dominated by a massless scalar, 
where (pit) ~ ^J^Mpi \nt, and a(t) ~ t^^^. This ends when the kinetic energy of the 
scalar redshifts down to the energy density of the black hole gas, roughly at a time ti ~ 
{Mpi/(pQ){(pl/ Pbh oY^^- Since pbh o is below the Planckian density, wiU have rolled sev- 
eral units of Planck mass until black holes begin to trap it, but not a huge amount. The 
general dynamics with the black hole terms is not tractable analytically, however we can 
follow the approach towards the attractor at 0* as long as — 0=,, ^ Mpi. In this case we 
can approximate the cosh and sinh functions by exponentials, and find a solution 

0(t) 
ait) 

The dimensionful parameter /i is determined by the scalar speed at the instant when 
the black holes begin to take over from the scalar kinetic energy. This limiting behavior is 
consistent with the numerical integration displayed in Fig. (2), describing a slow approach 
of towards the attractor 0=,,. 



0--M„ log (1 + /it) 
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After some time in this regime, the black hole density will have redshifted to near the 
scale of the potential V (0). If the initial value of g was close to unity, so that {p) / {q) ^ 
the modulus will remain in the region of ~ ~ 0{Mpi) where the minima of V{(f)) are 
located. The modulus will slide down the slope of V (0) gently, eventually settling down in 
the minimum and starting an inflationary era (if the value of V at the minimum is nonzero, 
or if the approach to the minimum is along a very fiat slope, where < H). We plot a 
typical evolution of the modulus towards the minimum in Fig. (3). This kind of behavior is 
generic whenever the initial black hole density is not too small, and when 0=,, ~ 0{Mpi), so 
that the black hole 0=,, point lies in the basin of attraction of a minimum of V. 




time time 



Figure 3: As in Fig. (2), but with nonvanishing potential V (cp) given by eq. (12). Again pbh o = 
10~^ Mpi. At the time ~ 10^ the modulus gets trapped in the local minimum of V located at 
0min ^ 5.8 Mpi. 

In the case when the initial value of pbh is too small, or if the attractor 0^, is far too 
far in the weak coupling regime, the modulus will not be trapped in the minimum of the 
potential V (0). Either the potential will become important too early, before the black hole 
gas manages to slow down the modulus, or 0* will lie outside of the basin of attraction of 
the minima of V. In each of these cases the modulus will overshoot the minimum at 0iiiin 
and will forever continue rolling down the slope V (0) towards the weak coupling. In such 
situations the Brustein-Steinhardt problem cannot be avoided and inflation will not start. 
The system will eventually converge to a tracking regime (analogous to the one observed for 
exponential potentials in quintessence models). A typical representative of this behavior is 
displayed in Fig. (4). 

However in those cases when the universe is created very close to the Planck scale, as 
would typically occur in for example the approach advocated in [17], the initial black hole 
density will be significant. This is simply because of equipartition and equivalence principle. 
Hence the moduli would be captured quickly in the black hole attractor 0*. It still remains 
that in the regions of the landscape where the initial value of the coupling was too small, (p^ 
would end up outside of the basin of attraction and the modulus would still run off to the 
weak coupling eventually. However, this behavior is endemic for the small initial coupling 
regardless of the subsequent dynamics. If the coupling had started off too small, the modulus 
would never have had a chance of ending in the minimum of V^(0) anyway. It would always 
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Figure 4: Here we take pbh o = 10 ^ Mp^. Now the number of black holes is not sufficient to 
stabilize the modulus and 4> overshoots the minimum at (j)min- 



fall out of the basin of attraction of the minima of V. Thus this initial regime is of very 
little interest to start with. One has to make some choices of the initial conditions, and the 
whole purpose of the resolutions of the overshoot problem here and in [19, 20, 21, 26, 27] 
is to demonstrate that the favorable initial region may be likely: in this case, it is a finite 
interval of the vevs of 0. In the landscape, at very high scales where the influence of the 
potential is negligible, the initial value of may be anything, and so there will be many 
regions where ~ 0{Mpi). In these regions our mechanism will work in a natural way, and 
they will inflate. The a posteriori probability to be in such regions will then be large because 
they will become exponentially big after inflation. 

Once inflation starts, the gas of black holes will dilute exponentially quickly. By the 
time inflation is over, assuming it produces about 65 e-folds or so, the number density of the 
black holes will drop by a factor of about 10^°. If inflation happens near the GUT scale, say 
V^^^ ~ lO^^GeV, the final number density of black holes at the end of inflation will be 

riBH final < 10"'°T^ ~ IQ-^^bV' . (16) 

Mpi 

After inflation, the universe will expand roughly for another factor of 10^° until it reaches 
the present epoch. The number density of these primordial black holes today will therefore 
be no more than about ubh /ma/(today) <, 10~^^^eV^, leaving us with at most one such 
black hole per more than 10^^ present horizon volumes. In other words, these primordial 
black holes will become just as efficiently diluted as the primordial monopoles in the original 
formulation of inflation [1]. They will quietly go away after completing their task. 

To summarize, in this paper we have proposed a new method for resolving the Brustein- 
Steinhardt moduli overshoot problem [16] in string cosmology. We find that a gas of primor- 
dial black holes in the early universe, where the mass of black holes depends on the modulus, 
can provide a transient attractor for the modulus. When such black holes are produced they 
will trap the modulus temporarily, and keep it within the basin of attraction of the minima 
of the nonperturbative modulus potential V. As time goes on, the black hole density will 
redshift away, placing the modulus gently on the potential slope. Following this the modulus 
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will slowly settle into the minimum, where it can drive inflation. A key ingredient of our 
mechanism is that the universe should be created near the Planck scale, so that the heavy 
states are produced initially with non- negligible number density. Such an approach has been 
proposed recently in [17], and our mechanism may be a natural ingredient for helping the 
modulus stop in the inflationary valleys in the landscape. After inflation has begun, the 
density of black holes redshifts away to exponentially small numbers, just like the density of 
heavy monopoles in the early models of inflation. These black holes are pushed outside of 
the current horizon. Wc stress that our mechanism relies mainly on the presence of heavy 
states in the theory whose mass depends on the modulus. It may be possible to realize a 
similar mechanism in other regimes, using different massive states. It would be interesting 
to consider such mechanisms in other cosmological models, where for example the initial 
density of massive states could be small, but the minima of the potential lie at very large 
vevs, such as in the theories with low scale uniflcation. 
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